Boundary effects on flux penetration in disordered superconductors 
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We investigate flux penetration in a disordered type II superconductor by molecular dynamics 
simulations of interacting vortices. We focus on the effect of different boundary conditions on the 
scaling laws for flux front propagation. The numerical results can be interpreted using a coarse 
grained description of the system in terms of a non-linear diffusion equation. We propose a phe- 
nomenological equation for the front position that captures the essential behavior of the system and 
recovers the scaling exponents. 
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I. INTRODUCTION 



In recent years, the discovery of high temperature su- 
perconductors has generated a renewed interest in un- 
derstanding the magnetization properties of type ff su- 
perconductors jj]. The magnetization process is usually 
described in terms of the Bean model as the mag- 
netic flux enters into the sample from the boundaries, 
quenched disorder is responsible for the formation of a 
constant flux gradient. The Bean model provides a phe- 
nomenological picture of average magnetization proper- 
ties, such as hysteresis and thermal relaxation Q, but 
does not describe local space-time fluctuations. Indeed, 
recent experiments showed that these fluctuations are 
not only common, but also large, spanning several length 
scales: the flux line dynamics is intermittent, taking place 
in avalanches 0, and flux fronts are fractal 

A widely used modeling strategy to describe the fluc- 
tuations in the magnetization process consists in molec- 
ular dynamics (MD) simulations of interacting vortices, 
pinned by quenched random impurities pj-pd[| . With this 
approach it has been possible to model flux profiles IB | , 
hysteresis 61, avalanches [|o|,0 and plastic flow 
One of the aims of these studies 0] is to link the macro- 
scopic behavior, as described for instance by generalized 
Bean models, to the microscopic vortex dynamics. 

Recently, we have shown that the flux penetration due 
to interacting vortices in a disordered superconductor can 
be described by a disordered non-linear diffusion equa- 
tion |12|. The equation can be obtained performing a 
coarse-graining of the microscopic equation of motion of 
the vortices. In the absence of pinning, the equation re- 
duces to the model introduced in Ref. |l3f| . This model 
has been solved analytically to provide expressions for 
the dynamics of the front for different boundary condi- 
tions fl3|,|l4| . When quenched disorder is included in the 
diffusion equations, flux fronts are pinned in agreement 
with MD simulations. Varying the parameters of the 
equation, we observe a crossover from flat to fractal flux 
fronts, consistent with experimental observations. The 
value of the fractal dimension suggests that the strong 



disorder limit is described by percolation |Q. In the 
weak disorder limit, we recover the analytical results de- 
rived in Refs. p3|JT^ ]. Using this description, we can thus 
directly link continuum theories, for which analytical so- 
lutions are possible, to the microscopic equations used in 
MD simulations Q. 

In this paper we systematically analyze the effect of dif- 
ferent boundary conditions on the propagation of the flux 
front. A similar study was presented in Refs. fll3||T4| ] in 
the framework of non-linear flux diffusion in the absence 
of disorder. We show here that the results are in agree- 
ment with MD simulations for various different bound- 
ary conditions. Next, we analyze the effect of disorder 
by varying the pinning strength in the MD system. The 
results are then interpreted theoretically by means of the 
non-linear diffusion equation |l2|-|i"4|] . Finally, we propose 
a phenomenological equation for the front position that 
is able to capture in a simple way the behavior of the 
system, recovering the numerical results for the different 
boundary conditions. 



II. VORTEX DYNAMICS MODEL 

In an infinitely long cylinder, flux lines can be repre- 
sented as a set of interacting particles performing an over- 
damped motion in a random pinning landscape [p|-pd]|. 
The equation of motion for each flux line i is given by 



(1) 



where the effective viscosity is obtained from material pa- 
rameters as r = $0^2/ PnC 2 ■ Here, $0 is the magnetic 
quantum flux, c is the speed of light, p n is the resistivity 
of the normal phase and H C 2 is the upper critical field. 
The first term on the right hand side of Eq. ([!]) accounts 
for the vortex-vortex interaction and it is given by 



j(f) = ^/(s^mfi/xy 



(2) 



where the function K\ is a Bessel function decaying ex- 
ponentially for |r| > A, and A is the London penetration 
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length Um. The interaction is cut off at a distance 6A 
to improve computational efficiency. The second term on 
the right hand side of Eq. (|l|) describes the interaction 
between pinning centers, modeled as localized traps, and 
flux lines. Here, G is the force due to a pinning cen- 
ter located at R p , I is the range of the wells (typically 
I « A), and p = l,...,N p (N p is the total number of 
pinning centers). For the pinning force, we use the fol- 
lowing expression: G(x) = — foS(\x\ — f) 2 , for |a?| < 1 
and zero otherwise. In the present simulations we re- 
strict ourselves to the case T = (see Ref. [|l(| for the 
implementation of thermal noise in MD simulations). 

As we discussed in the introduction, we intend to study 
the effect of different boundary conditions on the flux 
penetration. We start with an empty system and con- 
centrate all the vortices in a small strip at the boundary. 
Due to mutual repulsion, the vortices will be pushed in- 
side the material, forming a flux front. There are several 
ways to implement the boundary conditions, correspond- 
ing to different experimental situations. Here, we will 
consider the following boundary conditions [T^Jl^| : 

• (A) Constant total number of vortices. Experimen- 
tally this corresponds to an external control of the 
magnetic flux. 

• (B) Constant vortex concentration at the bound- 
ary. This case corresponds to an external control 
of the magnetic field. 

• (C) Total vortex number increasing at constant 
rate. This represents an external control of the flux 
rate. 

• (D) Boundary concentration increasing at constant 
rate, corresponding to a constant field rate. 

As a word of caution, one should notice that boundary 
conditions can be more complicated in reality, due to 
complex surface barriers that oppose flux penetration. 
These are not considered here and the only surface bar- 
rier is provided by already entered flux lines. 



III. MD SIMULATIONS 

We perform MD simulations based on Eq. ([!]) and an- 
alyze the flux front propagation for different values of the 
pinning strength / . We typically use up to N p = 800 000 
Poisson distributed pinning centers of width I = A/2 in 
a system of size (L x = 800X,L y = I00A), corresponding 
to a density of no = 10/ A 2 . The number N of flux lines 
depends essentially on the boundary condition adopted 
in the simulation. The injection of magnetic flux into 
the sample is implemented as in Ref. p2| , concentrating 
at the beginning of the simulation all the flux lines in a 
small strip L' <C A, parallel to the y direction, and im- 
posing periodic boundary conditions in both directions. 



The front position is taken as the x coordinate of the 
most advanced particle in the system at different times. 

The case (A) , corresponding to a constant vortex num- 
ber, was studied in detail in Ref. |]l2[| , where we showed 
that the front position x p grows initially with time as 
^1/3 f or sma n times. Eventually the front position slows 
down and saturates to a value £ p which increases as the 
strength of the pinning centers fo is decreased. In par- 
ticular, the front pinning length £ p was found to scale as 

/ 1//2 jjjl]. Here, we compare this behavior with that of 
cases (B-D). 

We first consider the case fo = 0, corresponding to a 
clean superconductor (i.e., without defects) in order to 
clearly identify the front penetration law in the initial 
regime. In Fig. [j] we show that the front advances as a 
power law with an exponent that depends on the bound- 
ary condition. In particular, we find x p ~ i 1 / 2 for case 
(B), x p ~ i 2 / 3 for case (C) and x p ~ t for case (D). 
The presence of disorder can affect all these behaviors 
in a different way, depending on the imposed boundary 
condition. The case (B) is quite similar to the case (A) 
studied in Ref. [Tig] : after an initial transient, the front 
gets pinned and the pinning length scales with the pin- 
ning strength (see Fig. |^). In particular, the pinning 
length scales roughly as x p ~ l//o, as shown in the in- 
set of Fig. ^. For case (C), the front never gets pinned. 
The effect of disorder is only to slow down the dynamics 
(Fig. ||). A similar behavior is found for case (D). 

IV. NON-LINEAR DIFFUSION 

In Ref. (l^] , we have shown that the front propagation 
can be described by coarse graining the system and ob- 
taining a disordered non-linear diffusion equation. The 
equation that rules the evolution of the local vortex den- 
sity p(r, t) is 

= V(apVp-pF c ) + k B TV 2 p , (3) 

at 

where a = <E>§ / 4 and F c is a random friction force due 
to the pinning centers, with a typical value scaling as 
F c ~ fo^/no- 

For T = and f = 0, Eq. (|) can be solved exactly 
using scaling methods |13| , |l4|| . In particular the density 
profiles obey the equation 

p{x,y,t)=t- a h{xfiP), (4) 

where a and f3 depend on the boundary conditions and 
satisfy a + 2(3 = 1. For the cases considered : A) a = 1/3, 
(3 = 1/3; B) a = 0, /? = 1/2; C) a = -1/3, [3 = 2/3; and 
D)a = — 1, j3 = 1. These results are in perfect agree- 
ment with the numerical simulations reported in Fig. [l], 
since the exponent [3 describes the dynamics of the front 
position. 

The function h(u) also depends on the boundary con- 
dition and for the case (A) is given by h{u) = (1 — u 2 )/6 
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for u < 1 and vanishes for u > 1. The other cases are 
reported in Refs. Jl^,|lJ]. We check by MD numerical 
simulations that, in the presence of disorder, the density 
profiles are described by the non-linear diffusion equa- 
tion (^). In Fig. ^ we show that the profiles follow the 
/o = solution and then deform when pinning starts to 
dominate. 



V. FRONT DYNAMICS 

In order to understand in a simple way the effect of 
disorder on the front propagation for different boundary 
conditions, we can write an equation for the average posi- 
tion of the front. The approach is very similar in spirit to 
what is done for the imbibition of porous media [|l7[ ■ As 
discussed above, the front is driven by the density gradi- 
ent against the pinning landscape. The density gradient 
can be estimated simply as Vp ~ p(0,t)/x p , where x p 
is the front position and p(Q, t) is the boundary density. 
The typical pinning force can simply be taken as fay/no- 
Collecting these two contribution we write 

Tdxp/dt = ap(0 1 t)/x p - foy/n^ . (5) 

In order to close the problem we have to specify the be- 
havior of the boundary density, which clearly will depend 
on the particular boundary condition chosen. Let us con- 
sider the various cases: 

(A) When the total number of vortices is conserved, the 
boundary density decreases as the front advances. This 
can be explained by noting that the total number of vor- 
tices can be roughly estimated as M = p(0,t)xp, so that 
p(0,t) = M/x p . Inserting this into Eq. (^) we obtain 

d,Xp/dt = 1/Xp — g, (6) 

where <?a = fcu/no/a and time is expressed in units of 
T /{Ma). Eq. (|6|) admits an implicit solution as 

(g/M)t = arctanh( v /g/Ma;p)/v / .g/M - x p , (7) 

which can be expanded for x p <C 1/ ^Jg/M to give 
Xp ~ t 1 / 3 (short times), and for x p ~ 1/^/g/M (long 
times) yielding x p ~ y/M/g{l - 2 exp(-2(g/M) 3 / 2 i)), 
corresponding to a pinned front. This behavior is in 
agreement with the scaling found in numerical simula- 
tions (see Ref. @). 

(B) The case of a constant vortex density at the boundary 
is similar to case (A). The boundary condition is simply 
p(0,t) — po, and Eq. (@) becomes 

dxp/dt = po/xp - 5, (8) 

where time is now expressed in units of T /{p^a). As in 
the previous case, Eq. (||) can not be solved explicitly 
but from the implicit solution it is possible to obtain 
the asymptotic behavior: x p ~ t 1 ! 2 at short times and 



Xp ~ po/g{l — exp(— (g/po)t)) at long times, in agreement 
with the results presented in Fig. ^. 

(C) This case is similar to case (A), with the difference 
that the total number of vortices increases with time (i.e. 
M = ht). Due to this the front is never pinned. In ab- 
sence of pinning we recover the t 2 / 3 behavior observed in 
MD simulations. We can not find an analytical solution 
of the equation in this case and resort to numerical in- 
tegration. The results presented in Fig. || indicate that 
the front asymptotically grows as x p — C + At 2 / 3 , where 
A is the solution of the equation A 3 = 3/2(/i — gA 2 ). 
Thus for low pinning we expect that the coefficient A de- 
creases linearly with the pinning strength. This result is 
in agreement with the MD simulation (see Fig. ||). 

(D) As in case (C), the front is not pinned by the dis- 
order, which has the only effect of reducing the front 
velocity. One can compute the asymptotic velocity im- 
posing Xp — Vt and inserting this expression in Eq. (|^), 
that now reads 

dxp/dt = ht/xp — g, (9) 

where h is the rate of increase of the boundary density. 
Solving for V one obtains V — (\/(? 2 + ih— <?)/2. 

VI. DISCUSSION AND PERSPECTIVES 

In this paper we have analyzed the effect of different 
boundary conditions on the flux penetration in disor- 
dered type II superconductors. We have conducted a 
series of MD simulations of interacting vortices and in- 
terpreted the results in terms of a non-linear diffusion 
equation. In the limit of no disorder, the equation has 
been solved in Refs. [[l3],f[4| yielding solutions for the front 
propagation and the density profiles. Here we have shown 
that these results are in perfect agreement with MD sim- 
ulations. Moreover, we have found that the presence of 
pinning centers affects the behavior of the system and, 
depending on the boundary conditions, the front is either 
pinned or simply slowed down. To clarify these effects, 
we have introduced a simple equation of motion for the 
front position, in the same spirit of Washburn approach 
to imbibition |T7j |, Despite its simplicity, the equation 
captures the essential features of the front dynamics. 

An important question that still remains to be ad- 
dressed is the relevance of these results for experiments. 
We can associate each one of the boundary condition 
studied here to different experimental conditions, corre- 
sponding to the way the field is applied to the sample. In 
many cases, the presence of surface barriers for flux pen- 
etration could in principle modify the scaling behavior 
discussed here. We believe, however, that away from the 
surface this effect should not be important and the scal- 
ing would be recovered. In addition, thermal and quan- 
tum creep effects in general may lead to a slowly moving 
front, even when pinning is expected. It is possible to 
account for these effects in this framework considering 
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the diffusion term in Eq. (|3|) or adding a random noise 
term to the front propagation equation. In conclusion, 
we expect that our approach together with a systematic 
series of magneto-optical measurements performed un- 
der different magnetic field controls (e.g., see Refs. |^| 
will represent significant steps for the understanding of 
the flux penetration phenomenon in disordered supercon- 
ductors. 
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FIG. 1. The average position of the front plotted as a func- 
tion of time. The data has been obtained from MD simula- 
tions with different boundary conditions (for a definition see 
text) in a clean system (/o = 0). The curves increase as t a , 
where a depends on the particular boundary condition im- 
posed. 
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FIG. 2. The average position of the front, obtained from 
MD simulations for case (B), plotted as a function of time. 
The curves increase as t 1 / 2 and saturates at long times to a 
value depending on fo- From top to bottom, the curves corre- 
spond to f = 0.001, 0.002, 0.003, 0.004, 0.005, 0.006, 0.007, 
0.01, 0.015, and 0.02. In the inset, we show that the pinning 
length scales as f$ . 
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FIG. 3. The average position of the front for case (C) ob- 
tained from MD simulations. The solid line corresponds to a 
clean system and the dashed line to a disordered system. The 
front moves as x v ~ At' 2 ^ 3 , and A is reduced in the presence 
of disorder. 
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FIG. 4. The density profile at different times t obtained 
from MD simulations for boundary condition (A) in the pres- 
ence of disorder. At the beginning the profile follows the 
solution of the non-linear diffusion equation in the disorder 
free case and eventually deforms due to the action of pinning. 
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